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Abstract. We present an analysis of the X-ray light 
curves of the magnetic cataclysmic variable DP Leo us- 
ing recently performed XMM— Newton EPIC and archival 
ROSAT PSPC observations. We determine the eclipse 
length at X-ray wavelengths to be 235 ± 5 s, slightly longer 
than at ultra-violet wavelengths, where it lasts 225 s. The 
implied inclination and mass ratio for an assumed 0.6 Mq 
white dwarf are i = 79.7° and Q = M wd /M 2 = 6.7. We de- 
termine a new linear X-ray eclipse and orbital ephemeris 
which connects the more than 120000 binary cycles cov- 
ered since 1979. Over the last twenty years, the optical 
and X-ray bright phases display a continuous shift with 
respect to the eclipse center by ~ 2.1° yr . Over the last 
8.5 years the shift of the X-ray bright phase is ~ 2.5° yr" 1 . 
We interpret this as evidence of an asynchronously rotat- 
ing white dwarf although synchronization oscillations can- 
not be ruled out completely. If the observed phase shift 
continues, a fundamental rearrangement of the accretion 
geometry must occur on a time-scale of some ten years. 
DP Leo is marginally detected at eclipse phase. The upper 
limit eclipse flux is consistent with an origin on the late- 
type secondary, L x ^ 2.5 x 10 29 ergs s" 1 (0.20 - 7.55 keV), 



■^j- , at a distance of 400 pc. 

7 ~ ' , Key words: Stars: binaries: eclipsing - stars: cataclysmic 
variables - stars:individual: DP Leo - X-rays: stars 
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1. Introduction 

DP Leo is one of the strongly magnetic cataclysmic bina- 
ries of AM Her type, a so called polar. It was discovered as 
the first eclipsing polar some 20 years ago as optical coun- 
terpart of the EINSTEIN source E1114+182 (Biermann et 
al. 1985), and continuously observed with from the ground 
and the space with e.g. HST (Stockman et al. 1994) and 
ROSAT (Robinson & Cordova 1994, hereafter RC94). It 
was found to be a two-pole accretor based on the de- 
tection of cyclotron emission lines in field strengths of 
30.5 MG and 59 MG, respectively (Cropper & Wickramas- 
inghe 1993). A thorough timing study by Robinson & Cor- 
dova (1994) using ROSAT X-ray data combined with ear- 
lier optical data revealed evidence for an asynchronous 
rotation of the white dwarf in the system. Asynchronous 
polars form a very small subgroup of all polars. There are 
four out of currently known 65 systems which show a small 
degree of asynchronism of typically about 1% (Campbell 
& Schwope 1999). With an extra spin of the white dwarf in 
DP Leo of about 2°-2.5° per year, the degree of asynchro- 
nism is seemingly much smaller than in the other four ob- 
jects. However, the earlier results are based on a mixture 
of optical and X-ray data with not necessarily common 
origin on the white dwarf. 

DP Leo was chosen as Calibration/Performance Ver- 
ification target of XMM-Newton and was observed with 
all three X-ray telescopes in Nov. 2000. The spectrum de- 
rived from these observations was recently published by 
Ramsay et al. (2001). Using a multi-temperature model 
of the post-shock flow, they found evidence of a very mas- 
sive white dwarf in excess of 1 Mq . In order to address 
the question of asynchronism in DP Leo based on X-ray 
data alone, we performed a timing analysis of the new 
XMM-data in combination with archival ROSAT obser- 
vations (one published by Robinson & Cordova, a second 
one unpublished, Sect. ||). For proper measurement of the 
eclipse parameters we used segmented data, where indi- 
vidual segments were determin ed w ith a Bayesian change 
point detection method (Sect. 2.3). Our main results are 
presented in Sect. |§|, where the eclipse parameters, an up- 
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dated eclipse ephemeris and the accretion geometry are 
discussed. The question whether there is a positive detec- 
tion of the secondar y at X-ray wavelengths in the eclipse 
is discussed in Sect. |3.5| . 

2. Observations and Data reduction 

2.1. XMM-Newton EPIC 

DP Leo was observed using XMM-Newton on 22 of 
November 2000 for a net exposure time of 19949 s. DP 
Le o was detected in a ll thr ee EPIC detectors (Turner et 
al. |200l| , Striider et al. |200l| ). The thin filter was used and 
the CCDs were read out in full window mode. 

Before extracting source photons, the data were pro- 
cessed using the current release of the XMM-Newton Sci- 
ence Analysis System (version 5.1). Standard procedures 
of data screening (creation of an image and a background 
light curve) revealed time intervals with enhanced parti- 
cle background. These intervals were excluded from the 
subsequent timing analysis using an approach described 
below (See Sect. 2JJ). This reduces the accepted exposure 
time to 15827 sec with the EPIC-PN detector. The obser- 
vations were performed without any interruption, i.e. full 
phase-coverage of the P or b = 5388 s binary was achieved 
with an average exposure of ~150s per 0.01 phase unit. 

2.2. ROSAT PSPC 

The field of DP Leo was also observed with the ROSAT 
PSPC on May 30, 1992 (ROR 300169, PI: Cordova) and 
on May 30, 1993 (ROR 600263, PI: Petre) for a net expo- 
sure time of 8580 s and 23169 s, correspondingly. Results 
of the 1992 observations were presented by Robinson & 
Cordova (1994), the results of the much more extended 
observations of 1993 are unpublished. 

Although the net exposure time of the two ROSAT 
observations was larger than the binary period, in neither 
case was complete phase coverage achieved due to 
the close proximity of the periods of the satellite 
and of the binary. 

There are further X-ray observations reported by Bier- 
mann et al. (1985), and Schaaf et al. (1987), respectively, 
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Fig. 1. Phase-averaged X-ray light curves of the ROSAT 
and the XMM-Newton observations. Phase bins have a 
size of 0.01 phase units of the P or b = 5388s binary. 

with the EINSTEIN and EXOSAT satellites. We make use 
of the timing of the X-ray eclipses detected with Einstein, 
the EXOSAT data have a too low count-rate and are not 
used further in this paper. 

2.3. Timing analysis with a Bayesian change point 
detection method 

In order to study the abrupt changes of the X-ray count 
rate particularly at eclipse ingress and egress we per- 
formed a timing analysis of the datasets using a Bayesian 
change point detection methocf] developed by Scargle 



1998, 2000). This method is well suited for a statistical 
examination when the arrival times of individual X-ray 
photons are registered (see Hambaryan et al. 1999). It is 
superior to methods which work on binned data, since it 
requires no a priori knowledge of the relevant time-scale 
of the structure which will be investigated. 

The method is applicable to data that are known to 
originate from a nearly ideal Poisson process, i.e. a class of 
independent, identically distributed processes, having zero 
lengths of dead time. The data gathered in XMM-Newton 



1 In general terms, the change-point methodology deals with 
sets of sequentially ordered observations (as in time) and de- 
termines whether the fundamental mechanism generating the 
observations has changed during the t ime t he data have been 
gathered (see, e.g. Csorgo & Horvath, 1997). 



EPIC-PN and ROSAT PSPC observations allow the mea- 
surement of arrival times of individual X-ray photons with 
a resolution of 73.3 ms and 0.1 ms, respectively, a resolu- 
tion much smaller than the ingress and egress time scale 
which is of the order of seconds. The EPIC-MOS data 
cannot be used for the study of the eclipse length, since it 
provides a resolution of only 2.6 s. 



Scargle's (1998, 2000) method decomposes a given set 
of photon counting data into Bayesian blocks with piece- 
wise constant count-rate according to Poisson statistics. 
Bayesian blocks are built by a Cell Coalescence algorithm 



(Scargle 2000), which begins with a fine-grained segmen- 
tation. It uses a Voronoi tessellation^ of data points, where 
neighboring cells are merged if allowe d by t he correspond- 
ing marginal likelihoods (see Scargle [2000 ) . 

We repeat here the essential parts of the 
method, expanding upon particular modifications 
of the original method as used in the present appli- 
cation. Assume that during a continuous observa- 
tional interval of length T, consisting of m discrete 
moments in time (spacecraft's "clock tick"), a set 
of photon arrival times D (fcj, ti+i, U+ n ) ls regis- 
tered. Suppose now that we want to use these data 
to compare two competing hypotheses, The first 
hypothesis is that the data are generated from a 
constant rate Poisson process (model M{) and the 
second one from two-rate Poisson process (model 
M2). Evidently, model M\ is described by only one 
parameter 8 (the count rate) of the one rate Pois- 
son process while the model M2 is described by 
parameters #1, 62 and t. The parameter t is the 
time when the Poisson process switches from 9\ to 
6*2 during the total time T of observation, which 
thus is divided in intervals T\ and T%. 

By taking as a background information (/) the 
proposition that one of the models under consider- 
ation is true and by using Bayes' theorem we can 
calculate the posterior probability of each model 
by (the probability that Mk (k = 1, 2) is the correct 
model, see, e.g., Jaynes 1997[ ) 

Pr(D\M k ,I) 



Pr{M k \D,I) 



-Pr{M k \I) 



(1) 



Pr(D\I) 

where Pr(D\M k , I) is the (marginal) probability of 
the data assuming model M k , and Pr(Mk\I) is the 
prior probability of model Mk (k = 1,2). The term 
in the denominator is a normalization constant, 
and we may eliminate it by calculating the ratio 
of the posterior probabilities instead of the prob- 
abilities directly. Indeed, the extent to which the 
data support model M 2 over M\ is measured by the 



2 The Voronoi cell for a data point consists of all the space 
closer to that point than to any other data point. 
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ratio of their posterior probabilities and is called 
the posterior odds ratio 



_ Pr(M 2 \D,I) 
' lx ~ Pr{M x \DJ) 



Pr{D\M 2 ,I) 



Pr{M 2 \I) 



Pr{Mi\I) 



(2) 



_Pr(D\M 1 ,I)_ 
The first factor on the right-hand side of Eq. (2) 
is the ratio of the integrated or global likelihoods of 
the two models and is called the Bayes factor for 
M 2 against Mi, denoted by B 2 \. The global likeli- 
hood for each model can be evaluated by integrat- 
ing over nuisance parameters and the final result 
for discrete Poisson eve nts can be represented by 
(see, for details, Scargle 1998 , 2000 , Hambaryan et 
al. 1999) 



£>21 — 



1 



B(n + 1,771-71+ 1) 



— ^ B ( n i + 1, mi - n x + 1) 



xB(n 2 + 

where B is the beta function, 



(3) 

1,7712 ~n 2 + l)Ar , 
rij and m h (j = 1,2), 
respectively are the number of recorded photons 
and the number of "clock ticks" in the observation 
intervals of lengths T\ and T 2 . At is the time in- 
terval between successive photons, and the sum is 
over the photons' index. 

The second factor on the right-hand side of Eq. 
(^) is the prior odds ratio, which will often be equal 
to 1 (see below), representing the absence of an a 
priori preference for either model. 

It follows that the Bayes factor is equal to the 
posterior odds when the prior odds is equal to 1. 
When B 2 \ > 1, the data favor M 2 over Mi, and 
when £>2i < 1 the data favor Mi. 

If we have calculated the odds ratio 2 i, in favor 
of model M 2 over Mi, we can find the probability 
for model M 2 by inverting Eq. (|^), giving 

Pr(M 2 \D,I) = -^-. (4) 
1 + U 2 i 

Applying this approach to the observational 
data set, Scargle's (1998, 2000] ) method returns an 



array of rates, (9i, 9 2 , cp ), and a set of so called 
"change points" (ti, t 2 , r cp _i), giving the times 
when an abrupt change in the rate is determined, 
i.e. a significant variation. This is the most prob- 
able partitioning of the observational interval into 
blocks during which the photon arrival rate dis- 
played no statistically significant variations. 

We determined the timing accuracy of these change 
points through simulations. We generated 1000 data sets 
(photon arrival times) with one change point each. The 
data in the two segments obeyed Poisson statistics. Each 
simulated data set had approximately the same charac- 
teristics as the observed data in terms of number of reg- 
istered counts, spanned time, characteristic time scales of 
expected variations, and was analyzed exactly in the same 
way. The standard deviation of the distribution of change 
points was found to be At cp = ±(2 — 3) s, if phase- folded 
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Fig. 2. Determination of the eclipse length for the ROSAT 
1993 (upper two panels) and the XMM-Newton observa- 
tions (lower panel) of DP Leo. The upper panel shows 
the distribution of the reciprocals of the time intervals be- 
tween neighbouring photons, the two lower panels binned 
X-ray light curves (bin size 5 sec). Vertical lines indicate 
the change points in the Poisson process, the space be- 
tween them is our measured eclipse length. 



data are used. The uncertainty was larger, when data in 
original time sequence were used due to the smaller total 
number of photons involved. We adopted an uncertainty 
of 2.5 s for the observationally determined change points 
which were used to derive the eclipse length. 

In Fig. ||]we visualize the outcome of the process 
for the ROSAT observations performed in 1993 
and for the XMM-Newton observations (EPIC-PN 
data only). The top panel shows the distribution of 
the reciprocal of the time interval between neigh- 
bouring photons (1993 data), in the two lower pan- 
els the X-ray light curves binned in intervals of 5 
sec are plotted. The change points determined by 
our method are indicated by vertical lines. 

2.3.1. XMM-Newton EPIC PN 

As a first step, we applied the change point detection 
method to a background region free of any X-ray source. 



A.D. Schwope et al.: An XMM-Newton timing analysis of the eclipsing polar DP Leo 



5 



Table 1. Features of the X-ray light curve derived from ROSAT and XMM-Newton observations of DP Leo. For a 
given epoch and instrument we list the mean count rate in the phase interval 0.80-0.90, the length of the eclipse At oc i, 
the phase of the center of the bright phase interval 4>c , and the length of the bright phase A</>b- The center of the 
bright phase is interpreted as accretion spot longitude. 



Epoch 


Mission/Det 


CR 

Is" 1 ] 


Atecl 

[s] 


<t>c 


A0 B 


1992.4 
1993.4 
2000.9 


ROSAT/PSPC 
ROSAT/PSPC 
XMM-Newton/EPIC 


0.35 
0.50 
0.25 


237 ±5 
233 ±5 
237 ±5 


0.006 ± 0.006 
0.013 ± 0.006 
0.067 ± 0.006 


0.57 
0.57 



This allowed us to determine time intervals where the 
background showed no significant variation. These were 
regarded as good time intervals and further used for the 
timing analysis of DP Leo. 

We extracted ~ 3180 EPIC-PN photon events from the 
source, whose arrival times were corrected to the solar sys- 
tem barycenter using the "barycen" task, as implemented 
in SAS version 5.1. 

We used phase-folded data and data in original time 
sequence in order to determine different quantities. The 
length of the bright phase and the eclipse length were mea- 
sured in phase- folded data, the times of individual eclipses 
(for a period update) were measured in original time se- 
quence. 

The mean bright-phase count-rate, the eclipse length, 
bright phase center and length of bright phase are listed 
in Tab. ^, whereas the times of mid-eclipse of the individ- 
ual eclipses are listed in in Tab. § The times given there 
are barycentric Julian ephemeris days, i.e. they take into 
account the 14 leap-seconds introduced between the first 
and the last data point. Since leap-seconds were omitted 
by seemingly all authors in the past (all timings in the 
literature are given in HJD only), we computed for 
all eclipse times we found in the literature the leap-second 
correction and include those times in the table for consis- 
tency and future work. 

2.3.2. ROSAT PSPC 

We also performed a timing analysis of the ROSAT PSPC 
observations. In both cases we used a circular region with 
30" radius to extract source photons. In total 2705 and 
8385 source counts were extracted, respectively. The ra- 
dius chosen encompasses 85% of the events in the ROSAT 
point spread function. Photon arrival times were corrected 
to the solar barycenter, as implemented in the Extended 
Scientific Analysis Software System (Zimmerman et al. 
1998| ). 

These data were treated in the same manner as the 
XMM-Newton data and the corresponding results are 
also listed in Tabs. |l| and ||. 

3. Results and discussion 



3.1. X-ray light curves and mean spectra of DP Leo 

The phase-averaged X-ray light curves of the two ROSAT 
and the XMM-Newton observations (summed signal from 
all three cameras) are shown in Fig. 0. At all occasions 
the source showed a pronounced on/off behavior with the 
eclipse roughly centered on the X-ray bright phase. The 
eclipse was covered 3 times in 1992, 8 times in 1993, and 
2 times in the PN-observation (good time intervals only, 
one eclipse was excluded from the analysis due to high 
particle background). 

In the 1992 observation the source showed a pro- 
nounced flare at phase 0.2. The end of the bright phase 
was not covered, the length of the bright phase, however, 
was inferred by RC94 from contemporaneous optical pho- 
tometry. A pre-eclipse dip, likely due to the intervening 
accretion stream occurred centered at phase 0.94. Inter- 
estingly, this feature was never observed again, indicating 
a re-arrangement of the accretion geometry. 

The 1993 observation covered the X-ray bright phase 
completely (although marginally at the start) thus allow- 
ing to measure the length of the bright phase from X-ray 
data alone. The source displayed similar brightness dur- 
ing the two ROSAT observations. The eclipse appeared 
centered on the bright phase. 

In 2000, the shape of the X-ray bright phase appeared 
almost unchanged compared to the 1993 observation. The 
eclipse now was clearly off-centered with respect to the 
bright phase. The rise to the bright phase was somewhat 
less steep than the fall. Compared with the earlier ROSAT 
observations, DP Leo appeared fainter in the center of 
the bright phase. According to Ramsay et al. (2001) and 
Pandel et al. (2001) DP Leo was in a state of interme- 
diate accretion at the time of the XMM-Newton obser- 
vations, whereas it was in a high state at the time of the 
ROSAT observations. The comparison of published results 
combined with our own analysis shows that the situation 
might be different. 

For the PSPC observations of 1992, RC94 derive a 
bolomctric blackbody luminosity for an assumed distance 
of 260 pc of L bb ,boi = KnF bh = 1.4t^ x 10 31 ergs" 1 . Scal- 
ing to the more likely distance of 400 pc gives Lbb,boi = 
3.3 x 10 31 ergs _1 . RC94 used a geometry factor k = 2. 
Ramsay et al. (2001) used k = sec(i — (3) — sec(80° — 
100°) = 1.06 and a distance of 400 pc and derive 1.5 x 
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Fig. 3. Relation between mass ratio and inclination for an 
eclipse length of 225 s. 

10 31 ergs~ 1 with the EPIC MOS detectors, more than 
twice that value with the EPIC PN detector. Within the 
accuracy of the measurements and scaled to the same ge- 
ometry factors the luminosities of the soft components at 
both epochs agree with each other. 

Contrary to the PSPC observations in 1992, there is a 
clear detection of DP Leo above 0.5 keV in the PSPC ob- 
servation performed in 1993, which allows fitting of a two- 
component spectrum. With the spectral resolution pro- 
vided by the ROSAT PSPC, the spectrum is well reflected 
by a combination of a black-body and a bremsstrahlung 
component. We fixed the bremsstrahlung temperature at 
the typical temperature of fcTbr = 15keV. The bolomet- 
ric flux in the bremsstrahlung component thus derived 
was Fbr,93 = 2.9 x 10~ 13 ergcm~ 2 s _1 . Application of the 
same simple model to the EPIC PN data, and adding a 
Gaussian for the iron line at 6.7 keV, gives a fitted tem- 
perature of fclbr = 1 1 ± 6 keV and a bolometric flux of 
-Fbr.oo = 2.4 x 10~ 13 erg cm -2 s _1 , which again is not in 
contradiction to the former ROSAT measurements. We 
conclude that the X-ray observations do not indicate an 
obvious change of the mass accretion rate between the 
three epochs. 

3. 2. Timing of the X-ray eclipse 

Application of our method to the data of DP Leo al- 
lowed an accurate determination of the eclipse length at 
X-ray wavelengths. The measurements at all three epochs 
agree with each other within the claimed accuracy, rang- 
ing from 233 s to 237 s with a 5 second accuracy 



Table 2. Times of mid-eclipse of all eclipses measured in- 
cluding new ROSAT and XMM-Newton data. Individual 
times are leap-second corrected times at the solar system 
barycenter (BJED: barycentric Julian ephemeris day). 
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(i) x = 


X-ray; O = 
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(see Tab. |l|). The average eclipse length is somewhat 
longer than that deduced from HST/FOS observations. 
Figure 5 of the paper by Stockman et al. (1994) implies 
that the eclipse length at ultraviolet wavelengths is about 
225s (measured at half intensity). The difference in the 
length of the ultraviolet and X-ray eclipses is due to the 
fact, that the source of X-ray emission is closer to the sec- 
ondary star than the source of the ultraviolet radiation. 
The former originates from the hot accretion spot while 
the latter has contributions from the whole surface of the 
white dwarf. We regard the eclipse length determined in 
ultraviolet data as relevant for the mass determination. 
Using the relation between eclipse length (in phase units), 
inclination and mass ratio given by Chanan et al. (1976) 
and a eclipse length of 0.0209 phase units (225 s), the rela- 
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tion between mass ratio and inclination as shown in Fig. |3| 
results. 

At the given period, the mass-radius relation for late- 
type stars by Caillault & Patterson (1990) predicts a sec- 
ondary star with mass of only 0.09 M . Assuming a typ- 
ical white dwarf with M w( j = 0.6 Mg, the mass ratio is 
Q = M wc i/Af2 = 6.7 and the implied orbital inclination 
i = 80° (cf. Bailey et al. 1993). Ramsay et al. (2001) ar- 
gued for a high mass white dwarf near the Chandrasekhar 
limit on the basis of their spectral model applied to the 
XMM-Newton data. This would imply a slightly higher 
inclination of i > 82°. However, the comparison of 
white dwarf masses based on their model with dy- 
namically determined masses in the well-studied 
polars QQ Vul (M dyn = 0.54 M vs. Mxst = 1.30 M ; 
Catalan et al. 1999, Cropper et al. 1999) and AM 
Her (M dyn = 0.45 M vs. Af X fit = 0.74 M ; Schwarz et 
al. 2001, Ramsay et al. 2000) shows that the X-ray 
spectralmodel tends to predict a too high mass for the 
white dwarf. A high mass white dwarf seems unlikely to 
us given the good fit shown by Bailey et al. (1993) to their 
optical eclipse light curve. A massive white dwarf would 
have less than half the radius of the 0.7 M white dwarf 
used by Bailey et al. and would not give a comparably 
good fit. We use, therefore, as a baseline for our further 
analysis the standard value M w d = 0.6 M . 

The new determination of the eclipse length has a 
much higher accuracy than that by RC94, who give 
216 ± 18 s and demonstrates the benefit of the Bayesian 
change point method. RC94 derive an upper limit of ^22 s 
on the length of the eclipse ingress/egress phase. The 
binned eclipse light curves of Fig. ^ clearly show, 
that eclispe ingress and egress lasts much shorter 
than 22 sec in the observations perfomed in 1993 
and 2000 but the count rate is not sufficient to 
resolve ingress and egress. We therefore cannot 
derive strong constraints on the lateral extent of 
the X-ray emission region. For comparison, in UZ 
For and HU Aqr where the egress phases could 
be resolved by EUVE and ROSAT observations, 
respectively, these features last only about 1.3 s 
(Warren et al. 1995, Schwope et al. 2001a), cor- 
responding to a full opening angle of the X-ray 
emission region on the white dwarf of only 3°. 

We proceed by updating the eclipse ephemeris of DP 
Leo by using X-ray data alone. We disregard optical data 
for this purpose since it is shown for other polars (e.g. in 
HU Aqr, Schwope et al. 2001b), that the optical and X- 
ray emitting regions might be disjunct. This could result 
in a shift of the X-ray with respect to the optical eclipse 
by several seconds and would corrupt the period determi- 
nation. A linear regression to the eclipse times measured 
with EINSTEIN, ROSAT, and XMM-Newton yields 

BJEDx,oci = To iCC i + E x P cc \ 

BJED x ,cci = 2448773.21503(2) + E x 0.0623628471(5) 



for the barycentrically and leap-second corrected time of 
the X-ray eclipse center. This is not to be mixed up with 
the orbital period of the binary system, since we are mea- 
suring the eclipse of a small structure on the white dwarf 
surface which obviously is not fixed in the binary system. 
Our determination of P ec \ is consistent with that of RC94 
only at the 2.5cr level with our period being longer. One 
reason for the slight inconsistency could be the omission 
of leap seconds by RC94, another the use of different types 
of input data, X-ray and optical data by RC94, X-ray data 
alone by us. 

3.3. Spot geometry and true phase zero 

The derivation of the true binary period of the system 
needs independent information about the eclipse of the 
white dwarf center (not the spot on it!). There is one di- 
rect measurement in the literature (Stockman et al. 1994) 
available to us based on HST/FOS measurements. One 
can correct from the observed X-ray mid-eclipse times 
to the times of mid-eclipse of the white dwarf. We did 
this for the eclipse data which entered the determina- 
tion of P cc \ above. This correction is based on the fol- 
lowing parameters: i — 79.65°, Q — 6.7, M w d = 0.6, mass- 
radius relations for the white dwarf and the secondary by 
Nauenberg (1972) and Caillault & Patterson (1990), spot 
latitude 100°, spot longitude at the different epochs as 
listed in Tab. [|, the longitude at the time of the EIN- 
STEIN observation was -22°. Usage of these parameters 
gives the correct eclipse length and length of the bright 
phase, if a height of the emission region of 0.02 R w( j is 
taken into account. The assumed height is in accord with 
other polars (e.g. Schwope et al. 2001a). The small cor- 
rections to the eclipse times as listed in the above table 
are +5.7, -0.5, -1.2, -6.2s for the EINSTEIN, ROSAT 
1992, ROSAT 1993, and the XMM-Newton observations, 
respectively. A linear regression to the corrected timings 
including the HST measurement gives the ephemeris of 
inferior conjunction of the secondary 
BJEDxUV.orb = T 0iOr b + E x P orb 

BJED XU v,orb = 2448773.21503(2) + E x 0.0623628460(6). 
The binary period P or b is slightly shorter than P oc i, as ex- 
pected. Both zero points, To j0r b and To iec i, agree with each 
other, since the spot was at about longitude zero at the 
epoch, when our cycle counting starts. The introduction 
of a quadratic term, as discussed by RC94, gives a bad fit 
to the data with an (O — C) time of about one minute 
at the epoch of the XMM-Newton observations and is 
therefore ruled out. 

3.4. Spot longitude variations 

The large shift of the bright phase with respect to the 
eclipse is interesting as such. RC94 noticed that since the 
early observations in 1980 the bright phase was continu- 
ously shifted from negative longitudes to about zero Ion- 
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Fig. 4. Accretion spot longitude variations as a function 
of time. The last two data points are from this work, the 
other points are from RC94. The dashed line is a linear fit 
to all data. 



gitude in 1992 and they deduced a yearly shift of the spot 
longitude of 2.05°. The data collected in Tab. [I] imply 
that this shift might become even larger, A% = 2.5° yr" 1 . 
Fig. |I] shows the synopsis of all the bright phase center 
(spot longitude) measurements. The phase shift is con- 
tinuous and monotonic over the last 20 years. The spot 
started at negative longitudes, i.e. in the half-sphere away 
from the ballistic accretion stream and now approaches 
the more typical (natural ?) location at about 30° (see 
Cropper (1988) for a compilation of accretion spot longi- 
tudes). 

Spot longitude variations can be caused by changes of 
the mass accretion rate, by synchronization oscillations or 
by an asynchronously rotating white dwarf. Accretion rate 
changes would not imply a monotonic phase shift, they 
would imply a positive spot longitude at high accretion 
rate and a smaller longitude at low accretion rate. Since 
the accretion rate most probably did not change consider- 
able between the ROSAT and the XMM-Newton obser- 
vations, spot longitude are difficult to explain this way. 

Synchronization oscillations are predicted to occur 
once a locked state between the white dwarf and the sec- 
ondary star is reached (Campbell 1989, King & White- 
hurst 1991). So far no measurement could be performed 
in order to test the theory, the relevant time-scales and 
the amplitudes of these oscillations. The predicted period 
of small oscillations about the locked state is P osc ~ 25 yr 
(Campbell & Schwope 1999), i.e. of the order of the time 
base covered meanwhile by the observations. There is no 
indication of a reversal of the spot longitude migration 
implied by an oscillation scenario. We therefore tend to 
favor the scenario of a dis-locked white dwarf and thus 
add DP Leo to the small sub-class of asynchronous po- 
lars with so far four members only (Campbell & Schwope 



1999). If our assignment is correct, DP Leo is different 
from the other systems in this sub-class showing a much 
smaller degree of asynchronism. RC94 already estimated 
the deviation (P or b — -frot)/-forb — 1CT 6 , whereas the ab- 
solute of this quantity in the other four is ~ 10~ 2 . We 
note that we cannot properly measure the spin period of 
the white dwarf in DP Leo, since the accretion spot is 
not fixed in the magnetic coordinate system of the white 
dwarf. Should the degree of asynchronism be of the order 
as derived here, a fundamental re-arrangement in terms 
of a pole-switch must occur sometimes in the not too far 
future. 

3.5. X-ray emission from the secondary star? 

We searched the XMM— Newton data for photons in the 
eclipse, which would be ascribed to the putative active 
secondary star. Omitting the first and last 10 seconds of 
the eclipse the total exposure time in eclipse investigated 
by us was 1436 s and included 6 eclipses (three cameras) 
in good time intervals. In the source-plus-background re- 
gion 29 photons were registered, while in the neighboring 
background region only 18 photons were registered. 

In order to estimate the likely count rate only from 
the source, we employed a Bayesian estimate using a 
method described by Loredo (1992), which is applicable 
to a dataset with low number of counts having a Poisson 
distribution. 

The most probable value of the count rate was 
0.0075 cts s _1 taken from the evaluated full Bayesian prob- 
ability distribution function. A Bayesian credible region (a 
'posterior bubble') is 0.0030-0.0120 cts s" 1 in a 68% (la) 
confidence interval. The 99.73% (3er) credible region gives 
a value 0.0000 — 0.0212 cts s _1 , consistent with zero. We 
regard our finding as uncertain marginal detection of the 
secondary in X-rays. 

With the nominal count rate the luminosity of the sec- 
ondary is 

L x = 2.5 x 10 29 (£>/400pc) 2 ergs s" 1 (0.20 - 7.55keV). 
For this estimate we used a count to flux conversion fac- 
tor ~ 1.6 xl0~ 12 ergs cm" 2 s _1 cts" 1 , adopted from a 
spectral study of one of the M5 type stars available in the 
XMM— Newton Lockman Hole data (see also Hasinger et 
al |2001[ ). 

This estimate is consistent with coronal emission of 



late type stars from the solar vicinity (Hunsch et al. 1999) 



4. Summary and conclusion 

We have analyzed XMM-Newton observations of the 
eclipsing polar DP Leo performed in the Calibra- 
tion/Performance Verification phase in November 2000 in 
parallel with former ROSAT-PSPC observations. The cen- 
ter of the bright phase indicates a accretion spot longitude 
of ^24° at the epoch of the XMM-Newton observations. 
Compared to former observations, it was shifted towards 



A.D. Schwope et al.: An XMM-Newton timing analysis of the eclipsing polar DP Leo 



9 



later phase, a continuation of a trend over the last 20 years. 
This finding is regarded as indicative of an asynchronous 
white dwarf, although synchronization oscillations around 
an equilibrium position cannot be ruled out. The differ- 
ence between the binary period, P or b> and the white dwarf 
rotation period, P rot , is small (P orb - P ro t)/Porb - 10~ 6 , 
4 orders of magnitude smaller than for any other of 
the presently known four asynchronous polars. We have 
derived accurate ephemerides for the center of the X- 
ray eclipse and inferior conjunction of the secondary in 
DP Leo, corrected for the leap seconds introduced be- 
tween 1979.9 and 2000.9 and reduced to the solar sys- 
tem barycenter. This can be used as reference for further 
studies of the evolution of the period of the binary and 
of the spin period of the white dwarf. Although we were 
able to correct the observed X-ray eclipse times to true 
binary phase zero through eclipse modeling, a direct de- 
termination of the conjunction is highly desirable. This 
would allow the proper measurement of a possible spin- up 
or spin-down of the white dwarf in DP Leo. The only way 
to achieve this in DP Leo is through high-speed photom- 
etry (preferentially in the ultra-violet), thus revealing the 
white dwarf. The determination of conjunction of the sec- 
ondary star be e.g. a spectral tracing seems to be impos- 
sible due to the faintness of this low-mass star. It was not 
yet detected spectroscopically, Bailey et al. (1993) photo- 
metrically derive an P-band magnitude of 2L m 8. 

It is interesting to note, that only the spot longitude 
displays large-scale shifts but not the latitude. A shift in 
latitude would result in a different length of the bright 
phase. No such effect is observed, the length of the bright 
phase is always ^0.57 phase units. Should an oscillation 
scenario be applicable to DP Leo, the non-observed lati- 
tudinal shift is an important extra datum for theory. The 
model of e.g. King & Whitehurst (1991) predicts a large- 
amplitude out-of-plane oscillation (60° or more). 

Our detection of the secondary in X-rays is marginal, 
the derived flux and luminosity are in agreement with that 
of single M-stars, which rotate much slower. However, a 
secure statement about the X-ray flux of the secondary 
star requires a much deeper exposure. 
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